Equazione di Schrödinger

Erwin Schrödinger

In meccanica quantistica, l'equazione di Schrödinger è un'equazione fondamentale che determina l'evoluzione temporale dello stato di un sistema, ad esempio di una particella, di un atomo o di una molecola.

Formulata da Erwin Schrödinger nel gennaio 1926 e subito pubblicata, è un'equazione differenziale alle derivate parziali, lineare, complessa e non relativistica che ha come incognita la funzione d'onda $\psi$, introdotta basandosi sull'ipotesi di de Broglie dell'onda di materia. Secondo l'interpretazione di Copenaghen, il modulo quadro della funzione d'onda è legato alla probabilità di trovare una particella in una determinata regione spaziale.

L'equazione di Schrödinger, fondante quella che verrà chiamata dall'autore meccanica ondulatoria, ha avuto un ruolo determinante nella storia della meccanica quantistica, ad esempio permettendo di comprendere perché soltanto alcuni valori discreti di energia sono ammessi per l'elettrone nell'atomo di idrogeno.

Storia e sviluppo

Agli inizi del Novecento Max Planck, per risolvere il problema del corpo nero, sollevato da John Rayleigh, gettò le basi della futura meccanica quantistica invalidando uno dei principi cardine della fisica classica, secondo il quale un dato sistema può assumere valori di energia che variano con continuità. Planck, al contrario, ipotizzò che l'energia della radiazione elettromagnetica non potesse assumere valori arbitrari, ma soltanto determinati valori discreti, multipli di una quantità elementare detta quanto, proporzionale alla frequenza propria di vibrazione dell'oscillatore. Applicando tale concetto, egli risolse il problema proponendo che gli oscillatori che generano la radiazione elettromagnetica possano oscillare solo con un’energia data dalla relazione:

$E=nh\nu$

dove $h$ è una costante detta costante di Planck, $\nu$ indica la frequenza della radiazione considerata e $n$ è un numero intero positivo.

Fu Einstein a dare significato fisico all'ipotesi di Planck, teorizzando che la radiazione elettromagnetica sia formata da pacchetti indivisibili di energia $E=h\nu$, nel 1926 chiamati fotoni. Nel 1922 Arthur Compton scoprì che un quanto di radiazione che urta un elettrone si comporta come una particella dotata di quantità di moto, pur avendo massa zero.

Nel 1924 Louis de Broglie ebbe un'intuizione: se la luce può comportarsi sia come onda che come corpuscolo, allora una particella, come ad esempio l’elettrone, potrebbe comportarsi anche come un'onda. Riprendendo concetti tratti dalla teoria della relatività ristretta, de Broglie ricavò la relazione:

$\lambda =h/p$

dove $p$ è la quantità di moto della particella considerata e $\lambda$ prende il nome di lunghezza d'onda di de Broglie.

Su tale base de Broglie riuscì a ricavare la quantizzazione del momento angolare nel modello atomico di Bohr, fino a quel momento introdotta ad hoc, e a dimostrare che le orbite atomiche erano effettivamente ellittiche-circolari, in linea con le previsioni di Bohr-Sommerfeld per l'atomo di idrogeno. A conferma della sua brillante idea, dal 1926 si susseguirono diversi esperimenti che confermarono l'ipotesi del dualismo onda-particella. Il più importante fu l'esperimento di Davisson e Germer del 1927 nel quale i due scienziati americani osservarono la diffrazione degli elettroni attraverso una serie di fenditure, da cui si ebbe la conferma che le particelle hanno la facoltà di comportarsi come onde.

Sulla scia di tali risultati Schrödinger andò alla ricerca di un'equazione che descrivesse il propagarsi della onda di materia di de Broglie. Rifacendosi alla fisica classica, dove il comportamento di un'onda e le sue caratteristiche sono espresse da una funzione continua, lo scienziato austriaco ottenne un'equazione differenziale le cui soluzioni, le funzioni d'onda, restituivano quei numeri quantici cruciali per la risoluzione della struttura atomica di un elemento. L'equazione di Schrödinger è inoltre in grado di descrivere l'evoluzione di una particella libera, rappresentando l'analogo della seconda legge di Newton $({\vec {F}}=m{\vec {a}})$ in meccanica classica, utilizzata per ricavare una predizione sulla traiettoria che un dato sistema seguirà date certe condizioni iniziali; va sottolineato tuttavia che in ambito quantistico il concetto classico di traiettoria perde significato. Sebbene questa teoria, chiamata dall'autore meccanica ondulatoria, possa sembrare un'interpretazione diversa dalla meccanica matriciale di Heisenberg, Born e Jordan, Schrödinger dimostrò che i due procedimenti matematici sono equivalenti.

Equazione

L'equazione di Schrödinger dipende dalle interazioni fra le varie componenti del sistema. Nel caso più generale l'equazione è scritta come:

$i\hbar {\frac {\partial \Psi }{\partial t}}(\mathbf {r} ,t)={\hat {H}}\Psi (\mathbf {r} ,t)$

dove:

• $i$ è l'unità immaginaria;
• $\mathbf {r} =(x,y,z)$ è un punto nello spazio tridimensionale;
• $\Psi (\mathbf {r} ,\,t)$ è la funzione d'onda, cioè l'ampiezza di probabilità per differenti configurazioni del sistema;
• $\hbar$ è la costante di Planck ridotta, cioè divisa per $2\pi$;
• ${\hat {H}}$ è l'operatore hamiltoniano, che ha la proprietà di essere hermitiano.

Il valore medio dell'operatore hamiltoniano sullo stato $\psi (x)$ rappresenta il valore atteso dell'energia:

$\langle {\hat {H}}\rangle _{\psi }=\int \psi ^{*}(x){\hat {H}}\psi (x)\,dx$

Poiché gli operatori posizione e impulso sono operatori hermitiani, l'equazione per sistemi scleronomi si può scrivere nel limite non relativistico:

$\left[i\hbar {\frac {\partial }{\partial t}}+{\frac {{\hbar }^{2}}{2m}}\nabla ^{2}-{\hat {U}}(\mathbf {r} )\right]\psi (\mathbf {r} ,t)=0$

dove intervengono gli operatori impulso $-i\hbar \nabla =\hbar \mathbf {k}$ e energia potenziale $U$.

L'ipotesi di de Broglie

I risultati ottenuti con l'ipotesi di de Broglie portarono allo sviluppo della meccanica quantistica intesa come meccanica ondulatoria. Con de Broglie, si associa a ogni particella un pacchetto d'onda del tipo:

$\psi (\mathbf {r} ,t)=\int d\mathbf {k} \,A(\mathbf {k} )e^{i(\mathbf {k} \cdot \mathbf {r} -\omega t)}$

che si propaga con velocità di gruppo:

$v_{g}={\frac {d\omega }{dk}}={\frac {dE}{dp}}={\frac {p}{m}}=c^{2}{\frac {k}{\omega }}$

dove $\omega$ è la frequenza angolare o pulsazione intesa come quella centrale del pacchetto d'onde, $\mathbf {k}$ è il vettore d'onda che identifica la direzione di propagazione del pacchetto, $E$ è l'energia associata alla particella e $p$ il suo impulso lineare.

Una volta associato il pacchetto d'onda alla particella, era necessario scoprire quale equazione fosse in grado di descrivere l'evoluzione del pacchetto d'onda compatibilmente alla meccanica quantistica e interpretarne le soluzioni. In tal senso applicando l'operatore di D'Alembert al pacchetto d'onde si ottiene:

$\nabla ^{2}\psi (\mathbf {r} ,t)-{\frac {1}{c^{2}}}{\frac {\partial ^{2}}{\partial t^{2}}}\psi (\mathbf {r} ,t)={\frac {m^{2}c^{2}}{\hbar ^{2}}}\psi (\mathbf {r} ,t)$

tenendo presente la relazione relativistica dell'energia :

$E^{2}=c^{2}\left(p^{2}+m^{2}c^{2}\right)$

L'equazione scritta sopra è l'equazione di Klein-Gordon in cui appare un termine a secondo membro che è un termine di sorgente della particella con lunghezza d'onda di Compton. Per particelle con massa nulla come i fotoni, l'equazione di Klein-Gordon è una normale equazione di D'Alembert che descrive la propagazione di un'onda elettromagnetica. Formalmente, tale equazione è ottenibile mediante le sostituzioni:

$E\rightarrow i\hbar {\frac {\partial }{\partial t}}$

$\mathbf {p} \rightarrow -i\hbar \nabla$

Limite non relativistico

L'equazione di Schrödinger nel limite non relativistico si può dedurre dall'equazione di Klein-Gordon. Considerando per l'energia la serie di Taylor al primo ordine:

$E=mc^{2}{\sqrt {1+{\frac {p^{2}}{m^{2}c^{2}}}}}\simeq mc^{2}+{\frac {p^{2}}{2m}}+\dots$

e definendo la frequenza angolare $\omega '$ non relativistica ottenuta da:

$\omega ={\frac {E}{\hbar }}={\frac {mc^{2}}{\hbar }}+{\frac {\hbar k^{2}}{2m}}={\frac {mc^{2}}{\hbar }}+\omega '$

si ottiene che usando l'espressione dell'energia non relativistica l'equazione di Klein-Gordon diventa l'equazione di Schrödinger:

$\left(i\hbar {\frac {\partial }{\partial t}}+{\frac {\hbar ^{2}}{2m}}\nabla ^{2}-U(\mathbf {r} )\right)\psi (\mathbf {r} ,t)=0$

Equazione di bilancio e di continuità

Dall'equazione di Schrödinger:

$i\hbar {\frac {\partial }{\partial t}}\psi (x,t)=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi (x,t)+{\hat {U}}(x)\psi (x,t)$

deriva un'equazione di bilancio. Per dimostrarlo si può prendere il complesso coniugato ad ambo i membri:

$-i\hbar {\frac {\partial }{\partial t}}\psi ^{*}(x,t)=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi ^{*}(x,t)+{\hat {U}}(x)\psi ^{*}(x,t)$

dove si è supposto che il potenziale $U(x)$ è reale. Moltiplicando l'equazione di Schrödinger per $\psi ^{*}(x,t)$

$i\hbar \psi ^{*}(x,t){\frac {\partial }{\partial t}}\psi (x,t)=-\psi ^{*}(x,t){\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi (x,t)+{\hat {U}}(x)\psi ^{*}(x,t)\psi (x,t)$

e per $\psi (x,t)$ la sua coniugata

$-i\hbar \psi (x,t){\frac {\partial }{\partial t}}\psi ^{*}(x,t)=-\psi (x,t){\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi ^{*}(x,t)+{\hat {U}}(x)\psi (x,t)\psi ^{*}(x,t)$

e sottraendo membro a membro si ottiene a sinistra dell'uguaglianza

$i\hbar \left[\psi ^{*}(x,t){\frac {\partial }{\partial t}}\psi (x,t)+\psi (x,t){\frac {\partial }{\partial t}}\psi ^{*}(x,t)\right]=i\hbar {\frac {\partial }{\partial t}}\left(\psi ^{*}(x,t)\psi (x,t)\right)$

e a destra invece

$-\psi ^{*}(x,t){\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi (x,t)+\psi (x,t){\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi ^{*}(x,t)=-{\frac {\hbar ^{2}}{2m}}\nabla \cdot \left(\psi ^{*}(x,t){\vec {\nabla }}\psi (x,t)-\psi (x,t){\vec {\nabla }}\psi ^{*}(x,t)\right)$

dato che il termine proporzionale al potenziale si cancella. Ma chiamando la grandezza scalare

$\psi ^{*}(x,t)\psi (x,t)=|\psi |^{2}=\rho$

e il vettore

$-{\frac {i\hbar }{2m}}\left(\psi ^{*}{\vec {\nabla }}\psi -\psi {\vec {\nabla }}\psi ^{*}\right)=\mathbf {J}$

dall'uguaglianza del membro di destra e di sinistra si deduce l'equazione di continuità:

${\frac {\partial \rho }{\partial t}}+\nabla \cdot \mathbf {J} =0$

Interpretazione della funzione d'onda

Si poneva a questo punto il problema del significato da attribuire alla funzione d'onda e più specificamente alla quantità $\rho =|\psi |^{2}=\psi ^{*}(x,t)\psi (x,t)$. Inizialmente si pensò d'interpretarla nel modo più "intuitivo", cioè come la densità di materia contenuta nel volume infinitesimo $dV=dx\,dy\,dz$. L'equazione di continuità avrebbe rappresentato così la conservazione della massa. Ma questa interpretazione risultò scorretta. Stessa sorte ebbe il tentativo di interpretazione come densità di carica. Successivamente Max Born mise in correlazione il concetto di funzione d'onda con la probabilità di rinvenire la particella in un punto qualsiasi dello spazio basandosi sull'analogia con la teoria ondulatoria della luce, per la quale il quadrato dell'ampiezza dell'onda luminosa in una regione rappresenta la sua intensità.

Si può esaminare il seguente esempio.

Dall'esperimento della doppia fenditura si evince che gli elettroni, anche se presi singolarmente, interferiscono con entrambe le fenditure come onde, per le quali si sommano le ampiezze. Nel caso di due onde sull'asse x:

$\psi (x)=\psi _{1}(x)+\psi _{2}(x)\$

Mutuando dall'espressione della intensità luminosa nel principio di sovrapposizione in ottica ondulatoria ($|{\vec {E_{1}}}+{\vec {E_{2}}}|^{2}$), esprimiamo la probabilità associata alle due onde nell'attraversare una o l'altra fenditura come:

$P=|\psi _{1}(x)+\psi _{2}(x)|^{2}=|\psi _{1}(x)|^{2}+|\psi _{2}(x)|^{2}+\psi _{1}^{*}(x)\psi _{2}(x)+\psi _{1}(x)\psi _{2}^{*}(x)$,

che è espressione evidentemente diversa da quella che sarebbe la classica composizione della probabilità nei due stati fondamentali in cui le due fenditure siano aperte separatamente: ($|\psi _{1}(x)|^{2}$ + $|\psi _{2}(x)|^{2}$).

Vi è infatti, per lo sviluppo del quadrato di un binomio, la comparsa dei due termini a destra, chiamati termini di interferenza, i quali determinano una gamma di valori, compreso lo zero, che giustificano le frange di distribuzione elettronica osservate sperimentalmente. Al concetto di interferenza ondulatoria viene quindi sostituito quello di "interferenza di probabilità".

Da tali considerazioni nacque l'interpretazione della funzione d'onda come ampiezza di probabilità e del suo modulo quadro come densità di probabilità. Ne deriva che l'equazione di continuità viene a esprimere la conservazione della probabilità.

Soluzioni dell'equazione di Schrödinger

Proprietà delle soluzioni

La classe di funzioni accettabili come soluzioni dell'equazione di Schrödinger è data dalle funzioni complesse definite su $\mathbb {R} ^{3}$ quadrato sommabili. Infatti deve essere:

$\int \,|\psi (\mathbf {r} ,t)|^{2}d\mathbf {r} =N$

cioè l'integrale deve convergere a un numero $N$ finito. Questa costante viene normalizzata, $N=1$, per compatibilità con il significato probabilistico della funzione d'onda: l'integrale esteso a tutto lo spazio del modulo quadro della funzione d'onda rappresenta la certezza, uguale a $1$, che la particella si trovi in un punto qualsiasi dello spazio stesso. Quindi le funzioni accettabili come soluzione sono le funzioni che appartengono a uno spazio lineare complesso chiamato spazio di Hilbert.

Dato che l'equazione di Schrödinger è lineare, date due possibili soluzioni $\Psi _{1}(\mathbf {r} ,t)$ e $\Psi _{2}(\mathbf {r} ,t)$, allora anche la funzione d'onda somma $\Psi _{1}(\mathbf {r} ,t)+\Psi _{2}(\mathbf {r} ,t)$ è a sua volta una soluzione dell'equazione. Questa linearità, espressa nel principio di sovrapposizione, ha numerose importanti conseguenze sulla dinamica e sulle proprietà dei corpi microscopici che obbediscono alle leggi della meccanica quantistica. Supponendo infatti di avere un sistema che ammetta una proprietà con due possibili valori, come ad esempio lo spin lungo una direzione, allora una possibile soluzione dell'equazione di Schrödinger prevede che il sistema possa essere in una combinazione lineare di due funzioni d'onda, una in cui ad esempio lo spin è up, l'altro in cui lo spin è down. In un certo senso, è come se una particella fosse in contemporanea in diversi stati e non in uno soltanto come viene invece sempre osservata. Questo problema fu affrontato da Schrödinger nel celebre paradosso del gatto.

Soluzioni per una particella libera

L'equazione di Schrödinger per una particella libera che si muove in uno spazio con una sola dimensione si riduce a

$i\hbar {\frac {\partial }{\partial t}}\psi (x,t)=-{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}}{\partial x^{2}}}\psi (x,t)$,

Seguendo la quantizzazione canonica, questa equazione può essere ottenuta semplicemente prendendo l'Hamiltoniana di una particella libera:

$H(p,x)={\frac {p^{2}}{2m}}$

e sostituendo all'impulso la derivata rispetto alla coordinata spaziale $x$, $p\rightarrow -i\hbar {\frac {\partial }{\partial x}}$, e considerando il quadrato $p^{2}$ come la doppia applicazione della derivata, cioè la derivata seconda.

Una possibile soluzione dell'equazione di Schrödinger per una particella libera è semplicemente la funzione

$\psi (x,t)=\exp {\left({\frac {i}{\hbar }}\left(px-{\frac {p^{2}}{2m}}t\right)\right)}$,

dove $p$ è al momento un semplice parametro libero. Infatti:

${\frac {\partial }{\partial t}}\psi (x,t)=-{\frac {i}{\hbar }}{\frac {p^{2}}{2m}}\psi (x,t)$

${\frac {\partial ^{2}}{\partial x^{2}}}\psi (x,t)=-{\frac {p^{2}}{\hbar ^{2}}}\psi (x,t)$,

da cui si ricava l'uguaglianza membro a membro:

$i\hbar {\frac {\partial }{\partial t}}\psi (x,t)={\frac {p^{2}}{2m}}\psi (x,t)=-{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}}{\partial x^{2}}}\psi (x,t)={\frac {p^{2}}{2m}}\psi (x,t)$.

La soluzione più generale può essere ottenuta considerando la sovrapposizione di diverse soluzioni con differenti parametri liberi $p$, cioè il pacchetto d'onda in una dimensione:

$\psi (x,t)={\frac {1}{\sqrt {2\pi \hbar }}}\int dp\,\phi (p)e^{i\left(px-{\frac {p^{2}}{2m}}t\right)/\hbar }$

dove il fattore prima dell'integrale è dovuto alla corretta normalizzazione e l'integrale sostituisce la somma per un parametro continuo come $p$.

Essendo un'equazione differenziale al primo ordine nel tempo, l'equazione di Schrödinger deve essere accompagnata dalla condizione iniziale che stabilisca la forma della funzione d'onda a un certo istante. Ad esempio al tempo $t=0$ si può imporre che la funzione d'onda sia:

$\psi (x,t=0)={\frac {1}{\sqrt {2\pi \hbar }}}\int dp\,\phi (p)e^{ipx/\hbar }$

in modo che la sua evoluzione nel tempo esista determinata per ogni istante $t$. Secondo l'interpretazione di Copenaghen della funzione d'onda la quantità:

$P(x,t)dx=|\psi (x,t)|^{2}dx\$

rappresenta la probabilità che la particella si trovi nell'intervallo $x,x+dx$, se la funzione d'onda è stata correttamente normalizzata:

$\int _{-\infty }^{\infty }|\psi (x,t)|^{2}dx=1$

in modo tale che la probabilità di trovare la particella in qualche punto dello spazio (in questo caso siamo su una retta perché stiamo prendendo in considerazione solo il caso unidimensionale, ma tutto ciò vale anche nel caso tridimensionale), deve essere 1 con certezza. Abbiamo inoltre stabilito che le funzioni accettabili come soluzioni dell'equazione di Schrödinger siano le funzioni definite in un campo vettoriale complesso e siano quadrato sommabili, cioè sia sempre verificata:

$\int _{-\infty }^{\infty }|\psi (x,0)|^{2}dx<\infty$

Il fatto che si tratti di un'equazione lineare implica che possiamo considerare la sovrapposizione:

$\psi (x,t)=c_{1}\psi _{1}(x,t)+c_{2}\psi _{2}(x,t)\$

(dove $c_{1},c_{2}\in \mathbb {C}$ rende valido il principio di sovrapposizione) che è anch'essa soluzione dell'equazione di Schrödinger. Un'altra caratteristica delle soluzioni dell'equazione di Schrödinger è che mentre il modulo quadro della funzione d'onda è importante perché rappresenta una probabilità, la fase dell'onda non ha rilevanza fisica.

Valori medi nelle rappresentazioni dell'impulso e della posizione

Poiché il significato della funzione d'onda è probabilistico si può parlare di valore medio di una grandezza fisica. Il valore medio della posizione (unidimensionale per semplicità) nella rappresentazione delle coordinate è dato:

$\langle {\hat {x}}\rangle =\int dx\,{\hat {x}}|\psi (x,t)|^{2}=\int dx\,\psi ^{*}(x,t){\hat {x}}\psi (x,t)$

e più in generale una qualsiasi funzione di $x$:

$\langle {\hat {f}}(x)\rangle =\int dx\,{\hat {f}}(x)|\psi (x,t)|^{2}=\int dx\,\psi ^{*}(x,t){\hat {f}}(x)\psi (x,t)$

Il valore medio dell'impulso è invece per analogia con il caso classico:

$\langle {\hat {p}}\rangle =m{\frac {d\langle {\hat {x}}\rangle }{dt}}=m{\frac {d}{dt}}\int dx\,\psi ^{*}(x,t){\hat {x}}\psi (x,t)={\frac {\hbar }{2i}}\int _{-\infty }^{+\infty }dx\left({\frac {\partial ^{2}\psi ^{*}}{\partial x^{2}}}{\hat {x}}\psi -\psi ^{*}{\hat {x}}{\frac {\partial ^{2}\psi }{\partial x^{2}}}\right)$

Risolvendo l'integrando che è uguale a:

${\frac {\partial }{\partial x}}\left({\frac {\partial \psi ^{*}}{\partial x}}{\hat {x}}\psi -\psi ^{*}x{\frac {\partial \psi }{\partial x}}-\psi ^{*}\psi \right)+2\psi ^{*}{\frac {\partial \psi }{\partial x}}$

si vede che:

$\langle {\hat {p}}\rangle =\int dx\,\psi ^{*}(x,t){\frac {\hbar }{i}}{\frac {\partial }{\partial x}}\psi (x,t)$

cioè

$p={\frac {\hbar }{i}}{\frac {\partial }{\partial x}}$

che è la definizione dell'impulso nella rappresentazione delle coordinate e più in generale qualsiasi funzione di $p$:

$\langle {\hat {f}}(p)\rangle =\int dx\,\psi ^{*}(x,t){\frac {\hbar }{i}}f\left({\frac {\hbar \partial }{i\partial x}}\right)\psi (x,t)$

Nella rappresentazione degli impulsi il valor medio dell'impulso è semplicemente:

$\langle {\hat {p}}\rangle =\int dp\,\phi ^{*}(p){\hat {p}}\phi (p)$

con il significato che $|\phi (p)|^{2}dp$ rappresenta la probabilità che la particella o il sistema abbia momento $p$ determinato nell'intervallo $p,p+dp$. Il valor medio di una qualsiasi funzione di $p$ in questa rappresentazione è dato:

$\langle {\hat {f}}(p)\rangle =\int dp\,\phi ^{*}(p,t){\hat {f}}(p)\phi (p,t)$

Vediamo infine che la posizione nello spazio degli impulsi è:

$\langle {\hat {x}}\rangle =\int dp\,\phi ^{*}(p,t)i\hbar {\frac {\partial }{\partial p}}\phi (p,t)$

cioè

${\hat {x}}=i\hbar {\frac {\partial }{\partial p}}$

che è la definizione della posizione nella rappresentazione dell'impulso e più in generale qualsiasi funzione di $x$:

$\langle {\hat {f}}(x)\rangle =\int dp\,\phi ^{*}(p,t)i\hbar f\left(i\hbar {\frac {\partial }{\partial p}}\right)\phi (p,t)$

Da notare la simmetria delle due rappresentazioni. Vediamo che la posizione $x$ e l'impulso $p$ devono avere valore medio reale, poiché sono grandezze fisiche devono essere direttamente misurabili, cioè sono osservabili del sistema, per cui deve valere:

$\langle {\hat {x}}\rangle =\langle {\hat {x}}\rangle ^{*}$

e

$\langle {\hat {p}}\rangle =\langle {\hat {p}}\rangle ^{*}$

cioè in meccanica quantistica le grandezze fisiche di interesse, le osservabili, tra le quali la posizione e l'impulso, sono operatori hermitiani e questo si può verificare direttamente.

Se calcoliamo il commutatore tra la posizione e l'impulso nell'asse x:

$[{\hat {p}}_{x},{\hat {x}}]\psi (x,t)={\hat {p}}_{x}{\hat {x}}\psi (x,t)-{\hat {x}}{\hat {p}}_{x}\psi (x,t)={\frac {\hbar }{i}}{\frac {\partial }{\partial x}}{\hat {x}}\psi (x,t)-{\hat {x}}{\frac {\hbar }{i}}{\frac {\partial }{\partial x}}\psi (x,t)={\frac {\hbar }{i}}\psi (x,t)$

cioè si ottiene una delle parentesi fondamentali della commutazione in meccanica quantistica:

$[{\hat {p}}_{x},{\hat {x}}]={\frac {\hbar }{i}}$

che significa che non si possono misurare simultaneamente posizione e impulso in una direzione con precisione. Mentre:

$[{\hat {p}}_{y},{\hat {x}}]=0$

$[{\hat {p}}_{z},{\hat {x}}]=0$

Queste parentesi sono le parentesi fondamentali della meccanica quantistica ed esprimono il principio di indeterminazione di Heisenberg.

Soluzione generale dell'equazione di Schrödinger (per potenziali non dipendenti dal tempo)

L'equazione di Schrödinger è l'equazione fondamentale della meccanica quantistica. La sua risoluzione in molti casi non è possibile, tuttavia esiste un metodo di risoluzione generale che va sotto il nome di metodo di Fourier per la soluzione di equazioni differenziali, il quale permette di ottenere importanti informazioni sul sistema ed è un metodo generale applicabile in molti problemi fisici di interesse.

Esplicitando l'operatore hamiltoniano dell'equazione di Schrödinger unidimensionale:

$i\hbar {\frac {\partial \psi (x,t)}{\partial t}}+{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}\psi (x,t)}{\partial x^{2}}}-{\hat {U}}(x)\psi (x,t)=0$

l'equazione può essere fattorizzata per variabili, cioè la soluzione può scriversi:

$\psi (x,t)=T(t)\cdot X(x)$

dove $T(t)$ è una funzione che contiene solo la variabile temporale e $X(x)$ contiene solo la variabile coordinata. Sostituendo la seconda relazione nella prima si ha:

$i\hbar X(x){\frac {dT(t)}{dt}}+{\frac {\hbar ^{2}}{2m}}T(t){\frac {d^{2}X(x)}{dx^{2}}}-{\hat {U}}(x)X(x)T(t)=0$

Separando le variabili ai due membri, cioè dividendo ambo i membri per $T(t)\cdot X(x)$, si ottiene che entrambi i membri devono essere uguali a una stessa costante che si denota con $E$:

$i\hbar {\frac {dT(t)/dt}{T(t)}}=\left(-{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}X(x)}{dx^{2}}}+{\hat {U}}(x)X(x)\right)\cdot {\frac {1}{X(x)}}=E$

Quindi si hanno due equazioni separate:

$i\hbar {\frac {dT(t)}{dt}}=E\cdot T(t)\qquad -{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}X(x)}{dx^{2}}}+{\hat {U}}(x)X(x)=E\cdot X(x)$

La prima equazione si risolve subito:

$T(t)=C\cdot e^{-iEt/\hbar }$

dove $C$ è una costante: questa dà la dipendenza dal tempo della funzione d'onda $\psi (x,t)$. La seconda equazione è detta equazione di Schrödinger indipendente dal tempo, e ha la forma di un'equazione agli autovalori per l'operatore hamiltoniano:

$\left[-{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}}{dx^{2}}}+{\hat {U}}(x)\right]X(x)=E\cdot X(x)\qquad {\hat {H}}X(x)=EX(x)$

essa rappresenta gli stati stazionari del sistema, cioè gli stati che non dipendono dal tempo. Questa equazione si risolve trovando lo spettro degli autovalori che può essere discreto o continuo o nel caso più generale sia discreto che continuo e quindi degli autovettori associati in modo che nel caso più generale la funzione d'onda si può scrivere:

$\psi (x,t)=\left(\sum C_{E}X(x)+\int C(E)X(x)\,dE\right)e^{-iEt/\hbar }$

dove $C_{E},C(E)$ sono coefficienti dipendenti da $E$, i primi nel caso discreto e i secondi nel caso continuo, mentre $E$ rappresenta l'energia del sistema. Quindi l'operatore hamiltoniano fornisce la dipendenza temporale della funzione d'onda e permette tramite la soluzione dell'equazione di Schrödinger di risolvere il problema agli autovalori per l'energia.

Autofunzioni dell'energia

Supponiamo inizialmente che la soluzione dell'equazione di Schrödinger indipendente dal tempo abbia uno spettro discreto di autovalori, allora la funzione d'onda $\psi (x)$ può essere sviluppata come:

$\psi (x)=\sum _{E}C_{E}u_{E}(x)\$

dove i coefficienti $C_{E}$ sono unicamente determinati, infatti:

$\sum _{E}u_{E'}^{*}(x)\psi (x)=\sum _{E}C_{E}u_{E'}^{*}(x)u_{E}(x)=\sum _{E}C_{E}\delta _{EE'}=C_{E'}$

dove $u_{E}(x)$ sono ancora le autofunzioni dell'energia. L'unica restrizione è che la funzione d'onda deve essere normalizzata e quadrato sommabile:

$\int _{-\infty }^{\infty }dx\,\psi ^{*}(x)\psi (x)=\int _{-\infty }^{\infty }dx\,|\psi (x)|^{2}=1$

e di conseguenza anche le autofunzioni devono essere normalizzate:

$\int dx\,u_{E'}^{*}(x)u_{E''}(x)=\delta _{E'E''}$

L'equazione di Schrödinger dipendente dal tempo si scrive:

$i\hbar {\frac {\partial }{\partial t}}\psi (x,t)={\hat {H}}\psi (x,t)$

Utilizzando lo sviluppo della funzione d'onda $\psi (x)$, la soluzione è:

$u_{E}(x,t)=u_{E}(x)e^{-iEt/\hbar }$

quindi la funzione d'onda completa è:

$\psi (x,t)=\sum _{E}C_{E}u_{E}(x)\cdot e^{-iEt/\hbar }$

Se lo spettro è continuo invece si può sviluppare la funzione d'onda in termini di autofunzioni di energia nel modo:

$\psi (x)=\int C(E)u_{E}(x)$

dove $u_{E}(x)$ sono ancora autofunzioni dell'energia, e queste devono essere normalizzate:

$\int u_{E'}^{*}(x)\cdot u_{E''}(x)=\delta _{E'-E''}$

dove interviene la funzione Delta di Dirac. I coefficienti $C(E)$ sono automaticamente determinati, infatti:

$\int dx\,u_{E'}^{*}(x)\psi (x)=\int C(E)u_{E'}^{*}(x)u_{E}(x)=\int C(E)\delta _{E-E'}=C(E')$

quindi la funzione d'onda completa è:

$\psi (x,t)=\int C(E)u_{E}(x)\cdot e^{-iEt/\hbar }$

L'equazione di Schrödinger unidimensionale ammette una proprietà speciale. Non esistono degenerazioni dei livelli di energia: a ogni autovalore $E_{n}$ discreto corrisponde un solo autostato.

Approssimazioni dell'equazione di Schrödinger

• Approssimazione di Born-Oppenheimer, nota anche come approssimazione adiabatica, è una tecnica usata in chimica quantistica e nella fisica della materia condensata al fine di disaccoppiare i moti di nuclei ed elettroni (cioè per separare le variabili corrispondenti al moto nucleare e le coordinate elettroniche nell'equazione di Schrödinger associata alla Hamiltoniana molecolare). Si basa sul fatto che le tipiche velocità elettroniche sono molto maggiori di quelle nucleari.
• Approssimazione orbitalica, tecnica usata in chimica quantistica. Ogni elettrone viene considerato singolarmente come appartenente a un atomo idrogenoide e la carica nucleare $Ze$, carica che viene utilizzata per calcolare il termine relativo all'energia potenziale da inserire nell'equazione di Schrödinger, viene corretta utilizzando la carica nucleare efficace $Z_{eff}$. Quindi la forma semplificata della funzione d'onda, utilizzata per descrivere un atomo polielettronico, diviene una funzione del tipo $\Psi =\Psi (n_{1})\;\Psi (n_{2})\;\Psi (n_{3})\ldots \Psi (n_{x})$.
• Approssimazione WKB in meccanica quantistica è un'approssimazione semiclassica nella quale si impone che la funzione d'onda sia una funzione esponenziale che varia lentamente e quindi viene espansa in potenze della costante di Planck.

Metodi risolutivi di approssimazione

• Metodo variazionale, rappresenta nella meccanica e chimica quantistica, un approccio utilizzato per trovare approssimazioni all'autostato di minore energia (stato fondamentale) e ad alcuni stati eccitati. Le basi di questo metodo si fondano sul principio variazionale.
• Teoria perturbativa, metodo di calcolo che si basa sull'introduzione, nell'hamiltoniana, di una perturbazione, ovvero un potenziale così piccolo da giustificare uno sviluppo in serie di potenze.
• Metodo di Hartree-Fock, metodo ab initio usato in chimica computazionale basato sull'approssimazione, valevole nel caso di fermioni, che $N$ funzioni d'onda associate al sistema siano approssimabili tramite un singolo determinante di Slater.

Evoluzione relativistica

L'equazione di Schrödinger, essendo basata sull'Hamiltoniana classica, non può rappresentare il comportamento quantistico di particelle con massa nulla, come i fotoni. Un importante sviluppo risultò pertanto fondamentale per introdurre anche il formalismo della relatività speciale, e portò alle due equazioni, di Klein Gordon e di Dirac, che rappresentano, rispettivamente, particelle a spin $0$ (dette anche particelle scalari) e particelle a spin ${\frac {1}{2}}$.

Bibliografia

• (DE) E. Schroedinger Annalen der Physik 79 p. 361; Annalen der Physik 79 p. 489; Annalen der Physik 81 p. 734 (1926).
• Traduzione in italiano degli articoli sopra, su ipparco.roma1.infn.it.
• Lev D. Landau e Evgenij M. Lifšic, Meccanica quantistica. Teoria non relativistica, Roma, Editori Riuniti, 2010, ISBN 978-88-64-73208-4.
• Peter Atkins e Julio De Paula, Chimica Fisica, 4ª ed., Bologna, Zanichelli, 2004, ISBN 88-08-09649-1.
• Donald A. McQuarrie, John D. Simon, Chimica Fisica. Un approccio molecolare., 1ª ed., Bologna, Zanichelli, luglio 2000, ISBN 88-08-17640-1.

Voci correlate

• Funzione d'onda
• Meccanica ondulatoria
• Operatore hamiltoniano
• Operatore posizione
• Operatore impulso
• Operatore di evoluzione temporale
• Autostato dell'energia
• Equazione di Schrödinger non lineare

Altri progetti

Altri progetti

• Wikiversità
• Wikimedia Commons

• Wikiversità contiene risorse sull'equazione di Schrödinger
• Wikimedia Commons contiene immagini o altri file sull'equazione di Schrödinger